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Abstract 

Duality transformation, which relates a high-temperature phase to 
a low-temperature one, is used exactly to determine the critical point 
for several models (2D Ising, Potts, Ashkin- Teller, 8-vertex), as the self 
dual condition. By changing boundary condition, numerically we can 
determine the self-dual(critical) point of the Ashkin- Teller (or Gaus- 
sian) model. This is the first explicit application of the duality to the 
numerical calculation, with the use of boundary conditions. 

In this short note, we propose a new method to determine the 2D Gaus- 
sian critical point of quantum spin chains. Although the finite-size scaling 
method is a powerful tool to determine the critical point, difficulty may 
occur for some cases. This difficulty comes from the structure of scaling 
operators. By changing the boundary condition, we have the other struc- 
ture of operators. Therefore, selecting boundary conditions, we can use the 
preferable structure to determine the critical point. The obtained results 
are summarized in eq.(]20|) and Fig. 2. 

As an effective theory of the ID quantum spin systems, the following 
sine-Gordon model (in Euclidean space-time) has been studied 

S = ^ J drdx [(9,(/.)2 + (5^.0)2] + J drdx cos ^20, (1) 

where a is the lattice constant. The dual field 9{t,x) (in the Gaussian 
language) is defined as 

drcP{T,x) = -d,{zK0), d,(P{T,x) = dr{iK9). (2) 
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We make the identification cf) = (j) + V^tt, 9 = 9 + \f2Tx. There exists the 
C/(l) symmetry for the field 9 but the second term of eq.(|^) violates the 
[/(I) symmetry for 0. For the free field theory, the scaling dimensions of the 
vertex operators exp(ibm\/20) and exp(ibim\/2i?^) are "n? /2K and Km? /2, 
where the integer variables n and m are electric and magnetic charges in the 
Coulomb gas picture. 

After the scaling transformation a — > e'^'a, we have the following renor- 
malization group equations 

dK-^ 1 2 dy K 

These are the famous recursion relations of Kosterlitz. Up to the first order 
of y, we find that y is an irrelevant field for i^T > 4 and relevant for K < A. 
There is a separatrix 'i2K~^ — Sln/C"^ — = which separates the infrared 
unstable region from the infrared stable region, and on this separatrix, the 
Berezinskii-Kosterlitz-Thouless transition occurs. The 2D Gaussian fixed 
line lies on y = 0. For K < A and y ^ 0, y flows to inflnity. For y > 0, (i;^) 
is renormalized to 7r/\/2 as y ^ +oo and for y < 0, {(/)) ^ as y ^ — oo. 

The infinite two dimensional plane can be mapped to a periodic strip of 
width L by the conformal mapping w = {L/2tt) log z (z = t + ix). In the 
rest of this note, we consider the boundary effect of this strip system. 

First let us consider the following 1-D Hamiltonian with the periodic 
boundary condition Q 

H = Ho + X dxOi, (3) 





where Hq is the fixed point Hamiltonian and Oi{= Oj) is a scaling operator 
whose scaling dimension is xi. According to Cardy||l|, the following finite 
size dependence of excitation energies up to the first order perturbation is 
obtained 

27r / /27r\^i^^ \ 

AEn = — {xn + 2^AC„i„ — + • • • , (4) 



L \ " \L J 

where L is the length of the system, Xn is the scaling dimension of the 
operator 0„. And C„i„ is the operator product expansion (OPE) coefficient 
of operators On and Oi as 

Oi{z,z)On{0,0) =Cnlni-] f^j On(0,0)+---, (5) 
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in which hi and hi are the conformal weights of Oi {xi = hi + hi). From 
eq.(^, we have the following RG equation 

dX 

{xi - 2)A. 



dlnL 



When xi < 2(relevant), the second order phase transition occurs at A = 
0, whereas xi > 2 (irrelevant), the second term in eq.(Q) is the finite size 
corrections of the excitation energies of the critical systems. Up to the first 
order perturbation theory, we find that at the point A = the scaled gap 
LAEn does not depend on the system size and the scaled gaps for several L 
cross linearly at A = 0. 

On the other hand, when the OPE coefficient Cnin becomes zero in 
some reason, the above argument is insufficient and we must consider the 
second order term of A in eq.(0). In this case, the scaled gap LAEn may 
have an extremum at the point A = 0. In practice, this is not a preferable 
thing, because the point of extremum is sensitive to finite size corrections 
of irrelevant operators such as L_2i-2l (x = 4). 

In the sine-Gordon model (||), we substitute the operator \/2 cos \/2</) 
for Oi. In this case, there is no operator On with a nonzero value of 
{Ol^{zi)Oi{z2)On{zs)) (which is related with the charge neutrality condi- 
tions in the Coulomb gas picture. Note that the operators e^^'^^^ are not 
allowed.) Thus the OPE coefficient in eq.(^ is zero. This indicates that 
we cannot expect the simple behavior of the finite size scaling method. In 
addition, for the irrelevant scaling field (x > 2), the system is in the massless 
phase, so the finite size scaling method does not work to determine the fixed 
points. 

If we put artificially half magnetic charges m = ±1/2 in the system, the 
OPE relations are 

Oi{z,z)Ot/,{0,0) = ^(^-j (^-j O^/,(0,0) + ..., 
Oi{z,z)O"y,{0,0) = -^[-) y 0^/2(0,0) + •••, 



where 



(6) 



Oi = \/2cos\/2(/>, 
0^/2 = V2cos-^0, (7) 
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and there are non-zero OPE coefficients in eq.(^). 

For a physical example with half-odd magnetic charges, Alcaraz, Barber, 
and Batchelorl^] considered the 5 = 1/2 XXZ spin chain using Bethe ansatz 

L 

H = — y^JS^ S^_^_i + SjSj^i + ASjSj_^_i], 
j 

with twisted boundary conditions Sf^_^_i ± iS\j^^ = e^**(5f it iS\), Sf^j^^ = 
5*^. When $ = 0, this model corresponds to the Gaussian model with 
K = vr/ arccos(A), — 1 < A < 1. According to their numerical results, the 
twisted boundary conditions change the electric and magnetic charges as 

n ^ n, 

m m -\ . (8) 

27r 

Hence when the twist angle ^ is vr, half odd integer magnetic charges appear. 
Recently Fukui and KawakamiO] studied this model analytically and their 
results are consistent with eq.(^. However, since their studies were based 
on the integrability, the off-critical behaviors were not treated. 

To see what happens when the boundary condition is changed in the 
Coulomb gas picture, we review the case of the following action 

I j-co j-L n ^ f°° 

S=—K dT dx{d^ef + ^K dTdJ{T,0). (9) 
27r J-oo Jo v27r J-oo 

Here we write the action with the field 9 which is dual to (f) and we assume 
the periodic boundary condition 9{t,L) = 9{t,0). If we transform the field 
9 as 9{t,x) 9{t,x) — ^x/\/^L, then we can eliminate the second term 
of eq.(^) with the additional constant term but the boundary 

condition is changed as 9{t,L) = 9{t,0) — <I>/\/2, which corresponds to the 
defect line along the imaginary time. When <I> = 2A^7r (A^ is an integer), 
this is the periodic boundary condition. 

After the dual transformation (|2|), the action@ is transformed as 

^ = 2^ y_oo Jo + j drdrHr, 0). (10) 
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This shows that there exist magnetic charges =F<I>/27r at r = ±00. Thus we 
obtain the ground state energy as|0] 



L ' ' ' ^ " 2 V27ry 
and the conformal anomaly number changes as 



>), (11) 



c{^) = 1-12xq{^) = 1-q(^^^ K. (12) 

We denote the state corresponding to the vertex operator Vn^m = 6*^"^+*^™' 
as |n,m). Since there exists a magnetic charge ^/2tt at t = —00, we find 
the change of this state as 

|n, m)<j> = |n, m + $/27r)$=0) (13) 

and because there exists a magnetic charge — <I>/27r at r = 00, the conjugate 
state is 

m| = o(ra, m + <I>/27r|. (14) 

Hence we obtain @ 

E„,„(*)-i^„(0) = |(|. + f (™ + |)'). (15) 

or 

En,mm - Em = T ( 2K + T""' + ^"^2^1 ■ ^^^^ 

From this equation, we find that the state |n, 0)$ corresponds to |n,<I>/27r)o 
which has the excitation energy Enfl{^) — Eq{(^) = Enfl{0) — Eo{0), and the 
momentum n^/L. 

Note that Dotsenko and FateevllEll considered the similar situation 



1 poo rL / <!>' \ 



(t„,0), (t„-oo), (17) 



in which the additional charge exists only at r = cx) but not at r = —00. 
The change of conformal anomaly number is the same as eg. (|l^) , if we set 
$ = 2$[Q]. But the structure of scaling operators is not same with the case 
of eq.(^. In their case, |n,m)$ = |n, m)o, and the conjugate state changes 



5 



as $(n, m| = o{n,m + $ /27r| (which is consistent to eq.(p^), but in (^) 
the conjugate relation does not change. This may be only the difference of 
picture. 

In the case of ^ = tt, we have half odd integer magnetic charges effec- 
tively. In this case, |0, — l)7r(= |0,— l/2)o) and |0)7r(= |0, l/2)o) are degen- 
erate for free field theory. Introducing the perturbation term of eq.(|l]) and 
using the first order perturbation theory, we obtain the hybridized states 

1^^1)^ = ^(10,-1). + 10).) (18) 

whose parity is even, and 

IV'2). = ^(|0,-1).-|0)^) (19) 

whose parity is odd. (Note that only when $ = and vr, parity is a good 
quantum number.) Using the OPE(^), we obtain the finite size dependence 
of energy up to the first order perturbation as 



E,W-E„(0) = _ --2.A- - 



(20) 



Thus we find that the energy eigenvalues of these states cross linearly at 
A = 0. 

In this stage we consider the symmetry of the states ([T8|), ([l9|). In 
the Ashkin- Teller language, the half magnetic charge operator \/2cos (l)/V2 
{\/2 sin (f) / \/2) corresponds to the operator P = fi^/i^ {P* = /i^(T^)[^]. The 
sine-Gordon model(|I]) is invariant under the transformation 

</'^<^+^, 0^0, andy^-y, (21) 

and the operators \/2cos</)/\/2 and -v/2sin0/\/2 are transformed as 

V2cos</)/\/2 -\/2sin(/>/\/2, 

V2sm(p/V2 \/2cos0/\/2, (22) 

thus at the point y = the system has the self-duality 0. 
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To verify the above things numerically, we study the following 5 = 1 
quantum spin chain, 



i7 = - S{-iy){S]S]^, + S]S]^, + A5|5|+i). (23) 
i=i 

The effective action of this model is described as eq.(|l]). The whole phase 
diagram was shown in ref.[§. The transition between the dimer and the 
Haldane phases is of the 2D Gaussian type. Using the Lanczos method, we 
calculate energy eigenvalues of finite systems (L = 8,10,12,14). Figure 1 
shows the scaled gap behavior of L = 10, 12, 14 systems with the periodic 
boundary condition for A = 0.5. We can see a minimum of the scaled 
gap. In Fig. 2, we show two low lying energies of the subspace J2 = 
with the boundary condition Sf^_^_^ = —Sf, S^_^_^ = —Sf, Sf^_^_-^ = Sf, 
which correspond to Ei{tt) and E2{tt). We see the expected behavior (20) 
for this twisted boundary condition. The obtained Gaussian fixed points 
agree with those obtained by the other method §. The conformal anomaly 
number is calculated as c = 0.998 for the periodic boundary condition and 
c{tt) = —3.185 for the $ = vr twisted boundary condition. In table 1, we 
show some extrapolated scaling dimensions. These numerical values are 
consistent with eqs.(pT|), (p^), (15), (20). With this method, we can also 



determine the Gaussian fixed line in the massless XY phase |10|] and apply 
to the S = 1 spin chains with the single ion anisotropy[jll], [T^ ]. 
Lastly we remark the case of the following sine-Gordon model, 



S = / drdx 

2'kK 



(dr4>f + (da;4>f] + I drdx COS VScj) 

J 27rQ^ 



with the same operator structure of (|^). In this case, the three point function 

(e^*^'^*^^i)\/2 cos \/8(/)(2:2)e^*^''^^^^'')o is not zero, so it is enough to consider 
the periodic boundary condition. The transformation corresponding to (|2l| ) 
is 

vr 

</'^0 + ;^, O^e, and y ^ -y. 
Nomura and Okamotop^ applied the crossing of excitations to determine 



the Gaussian fixed line, based on the RG analysis of Giamarchi and Schulz|14] 
The author thanks K. Nomura for illuminating discussions and critical 
reading. He also acknowledges K. Okamoto for useful discussions. The com- 
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Center, Institute for Solid State Physics, University of Tokyo. 
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scaling dimension 


xifl = 1/2K 


xo,i = K/2 


a^o,i/2(= xo{tt)) 




0.1786 


1.410 


0.3497 


K 


2.799 


2.819 


2.798 



Table 1: Scaling dimensions at the critical point A = 0.5, 5 = 0.2524. Here 
we extrapolated the corrections from the irrelevant field L2L2I (x = 4). 
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Figure captions 



Fig.l: The scaled gap behavior of L = 10(A), L = 12(n) and L = 14(0) 
systems with the periodic boundary condition for A = 0.5. 

Fig. 2: The low lying energies of the L = 14 system with the $ = tt boundary 
condition for A = 0.5. Parity even state(£?i(7r)) is denoted as solid line and 
parity odd state(£'2(7r)) is denoted as dashed line. The crossing point is the 
critical point and its estimated value is Sc = 0.2524. 
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